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ABSTRACT 

Motivated by recent models involving off-centre ignition of type la supernova explosions, we 
undertake three-dimensional time-dependent radiation transport simulations to investigate the 
range of bolometric light curve properties that could be observed from supernovae in which 
there is a lop-sided distribution of the products from nuclear burning. We consider both a 
grid of artificial toy models which illustrate the conceivable range of effects and a recent 
three-dimensional hydrodynamical explosion model. We find that observationally significant 
viewing angle effects are likely to arise in such supernovae and that these may have impor- 
tant ramifications for the interpretation of the observed diversity of type la supernova and 
the systematic uncertainties which relate to their use as standard candles in contemporary 
cosmology. 
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1 INTRODUCTION 



The most popular progenitor model for the average type la su- 
pernova (SN la) is a massive white dwarf, consisting of carbon 
and oxygen, which approaches the Chandrasekhar mass by an as 
yet unknown mechanism, presumably accretion from a compan- 
_ ion star, and is disrupted by a thermonuclear explosion (see, e.g., 
. iHillebrandt & Niemevej|200(]| for a review). This general picture is 
■ supported by the similarity of their photometric and spectroscopic 
properties which, after calibrating their peak luminosity by means 
of distance-independent propert ies, turns t hem into very good dis- 
tance indicators for cosmology jRiess et al. 1998: Perlmutter et al.l 
Il999h . While, however, the number of super novae observed at 
cosmological redshifts is steadily increasing jTonrv et al.l l2003l : 
iRiess et allliooi : lAstier et aUbood : Iciocchiatti et al.ll2006h , thus 
reducing the statistical errors of the cosmological parameters de- 
rived from them, a sound theoretical understanding of these objects 
-justifying in particular the calibration techniques applied for dis- 
tance measurements - is still lacking. 

Advances in numerical methods and the constant increase 
of computational power have allowed the development of three- 
dimensional (3D) simulations of the explosion phase of SNe la 
(Reinecke et al. 2002; Gamezo et al. 2003; Ropke & Hillebrandt 
l2005l ; lR6DkeetalJl2006l) . These facilitate a consistent treatment of 
instabilities and turbulence effects which play a key role in the ex- 
plosion mechanism. This level of sophistication allows the models 
to gain a high predictive power. At the same time, they provide 
the tools to address asymmetry effects in the explosion phase. The 
full 3D information fro m simulations has s o far only been used 
to derive nebular spectra jKozma et al ]l2005h . while synthetic light 
curves w ere obtained from sphe rically averaged results of such sim- 
ulations telinnikov et alj2006b . 



As the flame starts out in the sub-sonic deflagration burning 
mode, it is subject to buoyancy instabilities leading to complex 
morphologies of the ash region. While a turbulent cascade interacts 
with the flame on smaller scales, on the largest scales burning bub- 
bles float towards the star's surface evolving into a mushroom-like 
shape and partially merging with each other. The first full-star SN la 
simulation ( Ropke & Hillebrandt. 2005.) demonstrated that this pro- 
cess does not lead to global asymmetries by favouring low-order 
moments in the flow. Therefore, although complex in structure, the 
ash of the deflagration phase may be volume-filling on average if 
ignited in an initi al flame structure th at was uniformly distributed 
around the centre jRopke & Hillebrandt,2005,.) . The flame ignition, 
however, is an uncertain proc ess and it is possible that it proceeds 
in a very asymmetric manner teuhle n et al1l2006l ; iHoflich & SteinI 
12002). The consequence thereof may be an asymmetric flame evo- 
lution and the most extreme example is a flame ignited on only 
one side of the star. Floating rapidly towards the star's surface due 
to buoyancy, it may co nsume only material of one side of the star 
tein gale & D"ursill2007) . In some cases such models fail to release 
sufficient amounts o f nuclear energy to overcome the gravi tational 
binding of the star jCalder et al.ll2004l; IPlewa et alj|2004l) . Other 
cases, however, may explode the white dwarf and lead to very 
asymmetric compositions of the ejecta. Such a model will be dis- 
cussed in SectionHTl 



Another possible source of asymmetries in the ejecta is the 
propagation of a detonation wave. In delayed detonation models, 
such a wave is hypothesised to trigger after a period of burning 
in the deflagratio n phase iKhokhlov 1991). As recently discussed 
by iMazzali et al.l 120071) . this class of model is particularly ap- 
pealing since it provides a scenario which could be able to ex- 
plain most SNe la. One possibility for a deflagration-to-detonation 
transition (DDT) is the onset of the distributed burning regime 
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jNiemever & Wooslevlll997h . when turbulence first penetrates the 
internal flame structure. This happens at low fuel densities and con- 
sequently the DDT may take place at the outer edges of the de- 
flagration structure. The detonation propagating inwards from this 
spot may not catch up with the rapidly expanding structures on 
the far side before the fuel densities fall below the burning thresh- 
old; in this case, an asymmetric composition of the ejecta is ex- 
pected and correspo nding 3D simulations have been presented by 
iRopke & Niemev e? ( 2007). It is possible, however, that the DDT 
happens at multiple locations washing out the asymmetries in the 
ejecta. In an alternative scenari o, the "Grav itationally Confined 
Detonation" (GCD) suggested bv lPlewa et all (|2004), an asymmet- 
rically ignited flame fails to unbind the star and the ash erupts from 
the surface. Still gravitationally bound, it sweeps around the un- 
burned core of the white dwarf and collides on the far side. The 
resulting compression of fuel has been suggested to trigger a det- 
onation (Plewa et al. 2004) which potentially p roduces asymmetric 
ejecta compositions. However, recent results bv lRopke et al.l f2007h 
indicate that triggering a detonation in this scenario is possible only 
with special ignition setups in two-dimensional simulations. Three- 
dimensional models seem to disfavour this mechanism, because in 
those models greater energy release in the deflagration stage ex- 
pands the star such that a strong collision of the ashes is prevented. 
In some cases the energy release during deflagration is even suffi- 
cient to unbind the white dwarf (an example of one such model is 
discussed in Section|4]of this paper). 

The asymmetries produced in delayed detonation scenarios 
merit further investigation and will be the subject of a forthcoming 
study. Here, we are motivated to investigate and illustrate the pos- 
sible effects of chemical asymmetries on bolometric light curves. 
To this end, we will consider first a set of toy models which ex- 
plore the effect of an off-centre distribution of nuclear ash. We then 
extend the discussion to a real hydrodynamical explosion model - 
an asymmetrically ignited pure deflagration model that produced a 
weak explosion. 

The radiative transfer calculations required to obtain the light 
curves were performed using a 3D, fully time-dependent Monte 
Carlo code. We begin, in Section[2l by briefly describing the opera- 
tion of this code. Then, in Section[3] we discuss our set of artificial 
toy models. In Section|4] we present the results obtained w ith a re- 
cent hydrodynamical explosion model teopke et alj|2007' ) and in 
Section [5] we discuss the implications of this model. Our findings 
are summarised in Section|6l 



2 RADIATIVE TRANSFER CALCULATIONS 

Monte Carlo methods have been successfully applied to the mod- 
elling of radiation transport in supernovae for more than two 
decades. lAmbwani & Sutherland! |l988i), ''nd several subsequent 
studies, have used Monte Carlo simulations of 7-ray propagation 
to compute energy deposition rates and 7-ray spectra as functions 
of time for SN la. Subsequently, the Monte Carlo approach was ex- 
tended to follow not only the 7-rays but also the subsequent emis- 
sion and scattering of radiation in other spectral regions, thereby 
allowing bolom etric light cur ves to be obtained dCappellaro et al.l 
1 19971 ). Recently. ILucvI j2005l) has presented further generalisation 
of the methods and demonstrated that they can be readi ly used 
for three dimensional modelling. Following ILucvI ^20051) . these 
methods have been employed in a variety of contemporary multi- 
dimensional radiative transfer computat i ons of relevance to the 
study of SNe (see, e.g.. lKasen et al.ll2006l ; lMaeda et alj2006h . 



All the light curve calculations presented in this paper were 
performed usin^ the 3D Monte Carlo radiative transfer code de- 
scribed by Sim (2007), which is closely based on that presented 
by Lucv (2005). A brief description o f the operat ion of the code is 
given below but we refer the reader to Isiml ( l2007h for full details. 

The code assumes that the ejecta are in homologous expan- 
sion, an excellent approximation for the entire time span for which 
any significant radiative flux is able to propagate and escape. The 
supernova model is specified on a 3D Cartesian grid which expands 
smoothly along with the homologous flow. For the toy models (Sec- 
tion O a 100'^ grid is adopted while for the full explosion model 
(Section|4ll a 128^ grid is used. The model specifies the initial den- 
sity in each grid cell, the initial fractional mass of ^"^Ni and, when 
needed, the combined fractional masses of all iron group elements. 
The mass density at later times is readily deduced from the assump- 
tion of homologous expansion. 

Following Lucy (2005), the Monte Carlo quanta begin their 
lives as pellets of radioactive material. These pellets are placed in 
the ejecta in accordance with the initial distribution of ^"^Ni. As 
time passes in the simulation, the pellets decay as determined by 
the half-lives of the ^®Ni and daughter ^^Co nuclei. When a pellet 
decays, it becomes a 7-ray quantum, the propagation of which is 
then followed in detail. By means of Compton scattering or pho- 
toabsorption, the 7-ray quanta are able to deposit their ene rgy - 
the nu merical treatment of these processes is discussed by ILucvI 
( I2OO5I) . Whenever a 7-ray quantum is destroyed, it is assumed that 
its energy thermalises and is re-emitted as ultraviolet, optical or in- 
frared (UVOIR) photo ns; it is there after termed an r-packet in the 
nomenclature used by |LucvH2003) . 

Bolometric UVOIR light curves are obtained from the be- 
haviour of the r-packets. Currently, the code uses a simple grey- 
opacity treatment to follow the propagation and scattering of the 
r-packets (the specific forms of the opacity used for the models we 
present are discussed in Sections[3]and|4l. In principle, light curves 
can be obtained by angular binning of the r-packets as they emerge 
from the computational domain. However, greater efficiency is ob- 
tained by employing volume-based Monte Carlo estimators jLucvl 
1 19991) to deduce emissivities from the r-packet trajectories and then 
using these to perform a formal solution of the radiative transfer 
equation (details of the speci fic Monte C arlo estimators currently 
used in the code are given in ISinjl2007h . Thus, the end products 
of the calculation which are used in this paper are bolometric light 
curves computed for distant observers who are located on specific 
lines-of-sight and who are at rest relative to the centre of the super- 
nova model. 



3 TOY MODELS 

In this section, we undertake a preliminary investigation of the ef- 
fects of a lop-sided ''''Ni distribution on light curves using a set of 
simply-parameterized, artificial toy models. We begin by describ- 
ing these models and the numerical results obtained from them 
(Section 13.1b and then discuss the general implications of such 
models (Section[3]2]l. 

3.1 Description 

To explore the role of an off-centre distribution of ^^Ni, a grid of 
simplistic toy models has been constructed. In each model, a to- 
tal mass of 1.4 M0 is adopted. The mass density is assumed to 
be uniform and a maximum velocity of Umax = 10^ km s~^ is 
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Figure 1. Slices through the y-z plane showing the distiibution of ^®Ni adopted in the models. The dark grey areas indicate the regions in which there is no 
^^Ni while the light grey areas are those containing a fractional mass / of ^^Ni. The adopted values of / are indicated in the figure. Av/vaan specifies the 
fractional offset in velocity of the ^^Ni centre-of-mass from the true centre-of-mass of the SN. All the models have a total mass of 1.4 Mq and a total ^''Ni 
mass of 0.4 Mq. All the models are symmetric under rotation about the z-axis. 



chosen. This simple density distribution makes it convenient to ex- 
plore off-centre distributions of ^''Ni without necessitating changes 
in either the geometrical shape of the ^®Ni blob nor the underlying 
total mass distributions. In more sophisticated SN la models, the 
density distribution generally does extend to velociti es greater than 
our a dopted Umax (e.g. in the well-known W7 model. iNomoto et al.l 
Il984h - however, in these outer higher velocity regions, the density 
has typically fallen off by an order of magnitude or more compared 



to the inner parts which makes the outer regions relatively unimpor- 
tant in grey-bolo metric light curve calculati ons such as employed 
here. In particular. IPinto & Eastmai] ( l2000ah demonstrated that the 
light curve computed from a uniform-density spherical model with 
Umax — lO'' km s"^ agrees very well with that obtained from the 
W7 model (see their figure 2). 

Inside our toy SN models, the ^'^Ni is also given a spherical 
distribution but its centre-of-mass is offset in velocity relative to 



4 Sim et al. 



that of the SN; this offset (A«, which will generally be specified as 
a fraction of Umax) is one of the parameters of the models. Through- 
out the region which contains ^''Ni, a uniform mass-fraction / is 
adopted - this is the second parameter of the model. In all the mod- 
els the total mass of ^''Ni is fixed to 0.4 M0 ; therefore / determines 
the volume of the region which contains ^^Ni. In total, five models 
(A - E) are presented here; the structure of these models is illus- 
trated in Figure [T] The parameters which differentiate the models 
(Au/umax, f) are indicated in the figure and also tabulated in Ta- 
ble[T] In all the toy models, a uniform grey-absorption cross section 
of 0. 1 cm'^ g~^ is adopted for the treatment of "bolometric" UVOIR 
pho tons while the 7-ray transport and deposition is treated in detail 
(see lSinj2007l for a full description). 

Light curves have been computed for each of the five models 
(A ~ E) for three interesting viewing directions: 

(i) for an observer located at a large distance along the same 
direction in which the Ni is displaced relative to the centre of the 
SN; 

(ii) for an observer located diametrically opposite the observer 
in case (i); 

(iii) for an observer lying along a line-of-sight perpendicular to 
the direction of displacement. 

For each of these viewing angles in each of the models, the time 
(tp) and magnitude (Afp) of maximum light have been extracted; 
these values are reported in Table[T] 

The results in Table [T] indicate that an off-centre distribution 
of ^''Ni could introduce significant angular dependence in the light 
curve properties. In the models, the light curve peak is both ear- 
liest and brightest when viewed along the direction in which the 
^''Ni is displaced. This behaviour is as expected since the mean 
optical depth from the ''®Ni to the edge of the ejecta is lowest in 
this direction; therefore, more energy packets escape more quickly. 
Correspondingly, when viewed from the diametrically opposite di- 
rection, the peak magnitude occurs latest and is dimmest. 

It can be seen that the viewing-angle sensitivity of Afp is sig- 
nificantly affected by Ai; but is insensitive to the adopted value 
of /. As Aw is increased, the variation of A/p with angle in a 
given model increases from about 0.5 mag for the models with 
Au/wmax = 0.1 to more than 1.5 mag for the most extreme model 
considered (Au/umax = 0.3). 

To examine the angular behaviour of the light curve properties 
further, one model (model C) has been studied in greater detail. For 
this model, light curves have been extracted for an additional eight 
viewing directions. These directions have been selected to give a 
uniform grid of viewing angles in cos 9 (9 being the angle between 
the line of sight direction and the direction along which the ^''Ni is 
displaced). Figures |2] and [3] show the variation of Afp and tp with 
cos 9, respectively. 

Figure |2] shows that the peak magnitude follows a near-linear 
dependence on cos 9. The best-fit straight line (which is drawn in 
the figure) is given by: 

A/p = -18.67 - 0.52 cos 6* . (1) 

The near-linear variation with cos 9 is important since it means that 
the distribution of Mp one obtains by selecting random lines-of- 
sight is close to a top-hat function; thus the model predicts that if 
one were to observe it from a random viewing direction one would 
be equally likely to measure any value of Mp between about -18.2 
and -19.2 with no bias towards the median. 

A similarly simple relationship is apparent in Figure[3] Again, 
a linear fit is an adequate description; the best-fit straight line being 
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Figure 2. The peak bolometric magnitude (Mp) as a function of the viewing 
direction (8) for light curves computed with model C (shown as diamonds). 
The solid line is the linear fit to the results (see text). 
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Figure 3. The time at which peak magnitude occurs (tp) as a function of 
the viewing direction (8) for light curves computed with model C (shown 
as diamonds). The solid line is the linear fit to the results (see text). 



given by 

tp = 15.7 - 2.9 cos 9 (days) . (2) 

The absolute values for tp obtained from the toy models are 
generally somewhat low compared to those deduced by fitting tem- 
plates to observed light curves (e.g. typically ~ 19 days found in 
a rece nt study o f light curves from the Supernova Legacy Sur- 
vey by IConlev e t al. 2006). This is most likely a consequence of 
the simple treatment of the radiation transport which we have 
adopted (namely, the use of a grey, time-independent scattering 
cross-section per gram), but it may also be, in part, a consequence 
of the chosen structure of the toy models. However, a full investiga- 
tion of the sensitivity of tp to the assumptions made in the radiation 
transport and the construction of the models goes beyond what we 
require for this study; instead we focus only on the differential ef- 
fects introduced by the departures from spherical symmetry. The 
relative behaviour of tp and Mp is discussed below. 
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Table 1. Light cui've properties for the toy models (A - E). At)/ v^ax. specifies the velocity offset of the ^^Ni centre-of-mass relative to the maximum velocity in 
the ejecta. / gives the mass fraction of ^''Ni in the Ni region. Mp and tp are the magnitude and time at which the light curve reaches its peak. The superscripts 
t, i and 1. indicate properties as viewed by a distant observer lying along the direction in which the ^"Ni is displaced relative to the SN centre-of-mass, 
anti-parallel and perpendicular to this direction, respectively. AMp = Mp — is the range in AMp as a function of viewing angle for the model. 



Model Ai;/-!;n,ax / mJ k4 tj; AMp 

days days days 



A 


0.1 


0.45 


-19.00 


12.0 


-18.47 


15.2 


-18.71 


13.5 


0.53 


B 


0.1 


0.70 


-18.95 


13.9 


-18.43 


16.4 


-18.69 


15.6 


0.52 


C 


0.2 


0.70 


-19.21 


12.7 


-18.18 


18.1 


-18.67 


15.7 


1.03 


D 


0.2 


0.90 


-19.19 


13.5 


-18.14 


18.7 


-18.64 


16.1 


1.05 


E 


0.3 


0.90 


-19.43 


12.0 


-17.85 


20.3 


-18.59 


15.9 


1.58 



3.2 Discussion of the toy models and Amett's Rule 

The results obtained with the toy models described above can be 
used to gain some insight into how observable light curve proper- 
ties behave. In Figure [4] the values of A/p are plotted against tp 
for all the light curves summarised in Table[T] In addition, the grid 
of light curves computed to map out the angle dependence with 
Model C are also represented (eight additional points). We empha- 
sise that, since all the models considered adopt the same total mass 
of ^^Ni, the variation in Mp amongst the different points is en- 
tirely the result of the differences in the distribution of Ni with 
velocity. Since there is no compelling reason to suppose that all 
observed SN la synthesise the same total mass of ^''Ni - indeed, to 
the contrary, the total mass of ^^Ni is likely a dominant parameter 
in distinguishing SN la ~ our results (Figure|4]l cannot be regarded 
as predictions or a model for the directly observed distribution of 
SN la light curve parameters. Instead, they should be correctly in- 
terpreted in terms of a possible differential effect which is sensitive 
to the observer line-of-sight and which would appear in observed 
distributions superimposed upon other effects whose origins may 
lie in the underlying distributions of SN la explosion properties. 

As mentioned in Section [JTl the range of Alp obtained from 
the grid of models is very wide, spanning more than 1.5 mag. 
This illustrates that the effects of a lop-sided Ni distribution 
could, in principle, be very significant at the precision level set 
by contemporary observational data (e.g. in a sample of SNe for 
which bolometric peak magnitu des are available, as presented by 
IStritzinger & LeibundgutI llOOSl . the error estimates typically cor- 
respond to 10 - 15 per cent in flux). The scale of the effects is 
broadly comparable to those established for other models involv- 
ing departures from sphericity. For example. iHoflichI jl99lh inves- 
tigated some observational consequences of ellipsoidal SN ejecta 
and showed that for quite modest departures from sphericity (axis 
ratios in the range 0.9 to 1.2), one finds a viewing angle depen- 
dence of the peak brightnes s on the scale of tenth s of a magnitude 
(see also lHowell et alj200ll and lWang et alj2003l for motivation of 
such a model in the context of spectropolarimetric observations). 
As one would expect, even greater angular variations ~ comparable 
to those obtained here ~ can be obtained from ellipsoidal models 
with more extreme axis ratios (.Sim.2007 ). Similar scales of light 
cu rve angular depend ence (up to ~ 0.25 mag) were also found 
by iKasen et"al] ( |2004 ) in a stu dy of a SN model with a geometric 
"hole" in the ejecta and by Kas en & Plewal l l2006h in a very recent 
study of a 2D GCD model. 

We note that the angular variation is so large that one can 
likely rule out the more extreme models (models C, D and E) as 
being representative of typical SNe conditions - this follows from 
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Figure 4. The peak magnitude (Mp) versus the light curve rise time (tp) 
from the toy models. There are in total twenty-three hght curves repre- 
sented in the figure. The triangles represent the eleven light curves com- 
puted from Model C (as discussed in Section 13.11 these have been calcu- 
lated on a uniform grid in cos S-space) and the heavy dot-dashed line is the 
description of the Model C results obtained by combining the fits given in 
equations 1 and 2. The crosses represent the light curves computed from 
Model A while the diamond symbols show those from Models B, D and E; 
there are three viewing angles (parallel, perpendicular and anti-parallel to 
the direction of Ni-offset) per model (the numerical values are given in Ta- 
ble 1). The three curves - dashed, solid and dotted lines - show the expected 
Arnett relation obtained from equations 3 and 4 adopting a = 0.8, 1.0 and 
1.2, respectively. The sh aded grey-region indicates the r ange of values for 
Mp predicted following Stritzinger &Leibundgul ilOOj) for the ^GNi mass 
adopted in the toy models. 



the observed tightness of the correlation in the plot of bolomet- 
ric flux versus recessi on velocity for nearby SNe la (the "Hub- 
ble diagram"; see e.g. IStritzinger & Leibundgu3l2005l) . However, 
these models do show that it is possible to arrange extreme geome- 
tries that could produce substantial deviations from the mean and 
provid e one possible explanat ion for some apparently anomalous 
events ( iHillebrandt et alj2007h . 

Since both Mp and tp are simple functions of the viewing 
angle (equations 1 and 2), it follows that there is an equally simple 
relationship between these two quantities for a given model, as can 
be seen in Figure |4] For the models with moderate to large values 
of / (Models B - E), the points are all fairly close to following the 
same line in the Mp-tp plane. Thus, from these models alone one 
might hope that combined measurements of both Mp and tp could 
be used to extract the true luminosity despite the effect of viewing 
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angle. However, the points obtained with model A lie significantly 
below this line, showing that the relationship is not universal. 

A point of particular interest relates to the interpretation of Ar- 
nett's Rule tArn ett 1982) which is commonly used to estimate Ni 
masses in stud i es of SNe la light curves (e.g.lArnettetal.lll985l: 
BranchI Il992l: IStritzinge r & Leibundgu 3 llOOaT I Stritzinger et al.l 



2006al : 1stritzinger et al.l i2006b: Howell e t al.ll2006h . Arnett's Rule 



states that at maximum light, the luminosity roughly balances the 
instantaneous rate of energy generat ion in the SN. Adopting the 
convenient notation of lSranc 17( 119921) this relation is stated as 



aR{tp)Mm 



(3) 



where U is the peak luminosity, R(tp) is the rate of energy gener- 
ation at the time of maximum light (tp) and Mni is the total mass of 
^^Ni. The proportionality factor, a, is expected to be of order unity. 
R{t) is given by 



R{t) = [6.41 e^*/«'« + 1.34 e-*/"=''^]10*^ ergs s'^ M'^ 



(4) 



where t is measured in days. The numerical coefficients have 
been computed using the li fetimes and decay energies from 
lAmbwani & Sutherland' ( 1988). 

It is well-known from several previous studies that Arnett's 
Rule can usually be used to derive a reasonable estimate for Al^i 
from the light curve. However, it is not strictly obeyed, particularly 
if there is Ni p resent in the outermost r egions of the SN ejecta, 
as discussed bv lPinto & Eastmaj ( [200o3) . Since Arnett's Rule re- 
lates only global quantities, for any class of SN model in which 
a viewing angle dependence of the flux appears (including, but not 
limited to, the off-centre models discussed here), an attempt to mea- 
sure the nickel-mass via Arnett's Rule must introduce a systematic 
error which depends upon the direction of observation. In this sec- 
tion, we will use the toy models to investigate how this systematic 
would behave if real SNe (or at least a subset of them) were to 
harbour lop-sided distributions of ^^Ni. 

First we consider the case where a Ni mass is deduced by the 
application of Arnett's Rule (with a chosen value, or range of val- 
ues, for a) to direct measurements of the observed peak magnitude 
and time of maximum light. To elucidate this case, we over-plot 
in Figure |4] the curves obtained by combining equations 3 and 4 
for three different values of a (0.8, 1.0, 1.2). It comes as no sur- 
prise that Arnett's Rule with a fixed value of a does not describe 
the co-dependence of Mp and tp obtained from the models - the 
Rule is concerned with relating global properties of different SNe, 
not the detailed angular dependence within single objects. How- 
ever, the correlation of A/p and tp deriving from Arnett's Rule is in 
the same sense as that obtained from the models. This helps to sup- 
press the systematic error one would introduce with an unknown 
viewing angle; for the more moderate models (A and B), the sys- 
tematic error introduced would be only around 0.15 to 0.2 mag. In 
exceptional cases, however, a much larger systematic error of up to 
about 0.5 mag is possible for the adopted Ni mass. 

Secondly, we consider the case where reliable measurement 
of tp has not been possible so that an estimate of A/ni must be 
obtained from Mp alone . Adopting the relationship discussed by 
IStritzinger & LeibundgutI ( l2005j) (their equation 7 which is derived 
from Arnett's Rule and an empirically motivated assumption for 
the value of tp), leads to an expected value of Mp — —18.56 ± 
0.16 mag for A/wi = 0.4 Mq. This range of A-fp is indicated in Fig- 
ure|4]by the grey shaded region; t hus all the points which fall within 
this band are consistent with the IStritzinger & LeibundgutI | |2005|) 
relationship. For all the models considered, a significant fraction of 
the possible viewing directions lead to light curves which lie out- 



side the range. Thus, if the range of light curve properties produced 
by the models were representative of those from real SNe, using 
that relationship between Mp and A/ni in the analysis of a sample 
of observed light curves could lead one to infer a wider spread of 
nickel masses than required. 

The toy models have demonstrated that it is possible to con- 
struct very simplistic scenarios in which a lop-sided distribution 
of ^^Ni can introduce significant angular dependence in the emer- 
gent radiation. The scale of the variation is comparable to that 
introduced by other types of dep arture from spherical symmetry 
( lHoflich|[T99ll ; iKasen et alj|2004h and can be comparable to, or 
larger than, typical observational uncertainties. However, the toy 
models are very simplistic and are derived with no reference to 
realistic explosion physics. Therefore, they merely illustrate pos- 
sible effects and one must appeal to 3D explosion models to judge 
whether such effects are likely to have a part to play in reality. The 
remainder of this paper is concerned with the analysis of results 
from one such model. 



4 A 3D EXPLOSION MODEL 



4.1 The model 



The 3 D explosion simulation 3T2d200 described bv iRQpke et al.l 
( |2007|) has been used as the basis for a model to explore the ef- 
fects of an off-centre ^''Ni distribution in a realistic case. This sim- 
ulation followed the flame evolution when ignited in a lop-sided 
teardrop-like shape (upper left panel of Figure [5]l. Such an ignition 
configuration is motivated b y recen t studies of the pre-ignition con- 
vective burning phase dKuhlen et alj|2006l) . which suggest that the 
flow pattern is dominated by a dipole at this stage. A consequence 
could be a lop-sided ignition as adopted here, where the majority of 
ignition kernels are located on one side of the star but with a slight 
over-shooting across the centre. In this case, the flame propagates 
in both directions, subject to buoyancy instabilities and accelerated 
by turbulence. One side of the flame structure, however, remains 
dominant (see the top left panel in Figure[5). Although the ash bub- 
ble on this side of the star sta rts to sw eep around the core (similar 
to what has been described bv |Plewa"et a l. 2004), it does not collide 
on the far side (see the panels in the second row of Figure |5} be- 
cause the energy generated by nuclear burning is sufficient to grav- 
itationally unbind the star and expand it before a collision occurs. 
The lower left panel of Figure [5] shows the end-stage of this evo- 
lution; subsequently the expansion is approaching homology. The 
composition of the ejecta is illustrated in the lower right panel of 
Figure [5] showing a pronounced asymmetric bubble of iron-group 
elements. The simulation was carried out on a 512^ cells moving 
Cartesian computational grid (Ropke 2005). 

The model used here for radiation transport adopts a uniform 
Cartesian grid with 128"' cells. The total mass density and the mass 
fraction of Fe-group elements in each cell was obtained directly 
from the hydrodynamic simulation by combining sets of 4^ cells 
from the original 512^ grid. The cell properties were extracted from 
the latest time to which the hydrodynamics calculation was fol- 
lowed (~10s). For times beyond this point, including the entirety 
of the time for which the radiative transfer is followed, the ejecta 
are assumed to be in homologous expansion. 

Since the hydrodynamical results did not provide detailed in- 
formation on the composition of the ash material, it has been as- 
sumed that initially all the Fe-group mass was composed of ^^Ni; 
this gives a total of 0.448 Mq of ^''Ni. Since this investigation 
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Figure 5. Evolution of the 3D simulation 3T2d200. The initial flame configuration is shown in the upper left panel; the blue isosurface con'esponds to the 
zero level set associated with the thermonuclear flame. The following four panels (t = 0.2 s,t = 0.4 s, t = 1.0 s, and t = 3.0 s; left to right, top to bottom) 
illustrate the evolution of the model. Again, the blue isosurface con'esponds to the zero level set associated with the flame where burning is active and gives an 
approximate interface between fuel and ashes where the densities are below the burning threshold (lO'^ gcm^'^). The star is indicated in some of the figures 
by the volume-rendering of the density (red through yellow to white). The lower right image illustrates the composition of the ejecta (the front hemisphere is 
cut away) at the end of the simulation {t = 10.0 s). Unbumed material (carbon and oxygen) is shown in green and the white/yellow/orange volume-rendering 
corresponds to the density of the iron group elements. 
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is concerned with the viewing-direction dependence of the hght 
curve, this assumption will not significantly affect our conclusions 
- a lower mass fraction of ^''Ni in the ash would merely reduce the 
total luminosity; it would only affect the angle-dependence if the 
mass fraction were to vary significantly with position; whether or 
not this could be the case must await more detailed nucleosynthesis 
modelling in 3D. 



4.2 Treatment of opacity 

In the calculations using the 3T2d200 model, a slightly more so- 
phisticated treatment of UVOIR opacity is employed than the uni- 
form grey cross-section adopted for the toy models. Following 
jMazzali & Podsiadlowski' '2006Y an opacity which depends on 
composition is adopted, the particular form used here being 



-19.0 



(7 = 0.26 (O.gXpe-grp + 0.1) cm^ g~i 



(5) 



where Xpe grp is the mass-fraction of iron-group elements, which 
varies throughout the model. This simple parameterisation takes 
advantage of the basic compositional information available directly 
from the hydrodynamics calculation (namely the mass-fraction of 
heavy elements) to attempt to account for the dominance of the 

iron-group elements in providing opacity. 

^As has been discussed elsew here (e.g. see IPauldrach et al. 

19961 and IPinto & Eastmanll2000bl for in-depth discussion or lSiml 
20071 for comments in the context of the code used here), the 
UVOIR-opacity in SN ejecta is mostly dominated by spectral lines. 
Given this, the grey-approximation adopted here is not valid in de- 
tail - in particular, it cannot account for the process of photon re- 
distribution in frequency-space which allows photons to escape the 
ejecta by down-scattering into frequency regimes in which the den- 
sity of spectral lines is low. This process effectively limits the max- 
imum trapping any photon can undergo, an effect which is most rel- 
evant at the earliest times. For this reason, the grey-approximation 
may overestimate the dependence of radiative transfer effects on 
geometry. Nevertheless, the grey approximation is retained here 
since it dramatically reduces the computational demands of the ra- 
diation transfer calculations compared to a full non-LTE calcula- 
tion. Despite this approximation, the computations can be expected 
to correctly identify the sense in which purely geometrical effects 
will act and to give a simple estimate, or at least a reliable limit, for 
their scale. 

Of course, fully 3D, non-LTE, non-grey calculations are nec- 
essary to quantify the effects rigorously - however, to utilise this 
level of sophistication in the radiative transfer fully, one would 
also require more complete information on the detailed nucleosyn- 
thesis products and the ultimate evolution toward homology than 
is present ly ava i lable f or the 3T2d200 model. Fortunately, as dis- 
cussed bv lLucvl ( |2005|) . the Monte Carlo radiative transfer method 
is well-suited to incorporating more advanced treatments of micro- 
physics and thus there i s great promise for more detailed treatme nts 
when merited (see, e.g. lKasen et al.ll2006l . lKasen & Plewall2006[) . 



4.3 Results obtained with the explosion model 

We begin by presenting and describing the light curve properties 
obtained from the 3T2d200 model. Figure|6]shows light curves ob- 
tained for the special viewing directions parallel and anti-parallel 
to the vector connecting the centre-of-mass of the ejecta and the 
centre-of-mass of the ^®Ni; throughout the discussion we denote 
this vector n - note that n is nearly, but not quite, parallel to z. 
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Figur e 6. Computed light curves from the 3T2d200 model iRopke et alj 
1200711 . The dashed line shows the angle-averaged light curve. The heavily 
and lightly drawn histograms show, respectively, the light curve as seen by 
distant observers who lie along the direction-vectors n and —fi. fi is the 
unit- vector which points from the centre-of-mass of the ejecta to the centi'e- 
of-mass of ash material. 



We see from the figure that there are differences in the light 
curve properties as seen from these two directions - in agreement 
with the results from the toy models, the scale of the effect is sig- 
nificant compared to typical observational accuracy (^C 0.1 mag). 
Also in accordance with the results obtained in Section lJTl the light 
curve is brightest and rises to maximum fastest when viewed down 
n. In contrast to the toy-model results, however, the anti-parallel 
direction does not show the dimmest light curve for the 3D model 
- this is a consequence of the more complex distribution of ash 
material which determines both the location of energy generation 
(^^Ni) and the distribution of opacity (via the dependence on Xpe-gip 
in equation 5). As might be anticipated from Figure [5] (lower pan- 
els), the dimmest light curves are actually seen by observers whose 
line-of-sight is nearly perpendicular to n (see below). 

While the light curves shown in Figure |6] indicate that it is 
possible to see significant differences in the peak magnitude as a 
function of direction, it is important to consider how probable a 
randomly aligned observer is to see a given range of light curve 
parameters. To this end, a grid of 100 viewing directions (evenly 
spaced in solid angle), has been employed to map out the angular 
distribution of the light curve properties. In Figure|7]the peak mag- 
nitudes are shown for this grid of directions as a function of the 
polar angle, 9, (6 being the angle between the z-axis of the model 
and the observer's line-of-sight). There is clearly a net trend in the 
variation of Alp with 6; for large values of cos the mean trend is 
very similar to that found in the toy models. But, as noted above, 
owing to the complex distribution of ash material, a much weaker 
trend is revealed in the region where cos 6^ < 0. 

The grid of light curves as a function of viewing directions can 
be used to determine the probability distribution for the peak mag- 
nitude. This has been done and the resulting cumulative probability 
distribution is shown in Figure [8] 

Although more complex in detail, the probability distribution 
function (pdf) obtained shares most of the important properties sug- 
gested by the toy models; namely, it is fairly wide (spanning ~ 0.5 
mag) and the probability is not strongly concentrated around the 
median (for example, fully one quarter of the directions give peak 
magnitudes that are brighter than the median by at least 0.1 mag). 
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Figure 7. The peak magnitude as a function of tlie splieiical polar angle as 
obtained with the 3T2d200 model. For each value of 8 there are 10 points, 
each representing a different azimuthal angle, 
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Figure 8. The peak-magnitude cumulative probability distribution function 
- the probability that the obtained value of A/p is less than m - as obtained 
from the 3T2d200 model. The median value is -18.67 (dashed line). 

An important difference from the pdf implied by the toy models 
is a distinct asymmetry (as could be anticipated from Figure |7J - 
there is a probability tail extending to brighter magnitudes meaning 
that one may see larger differences on the bright side than on the 
dim side of the median. The probability distribution for A/p will 
be used in Section [F!2l to quantify the relevance of this work to the 
interpretation of the local Hubble diagram. 

For completeness, Figure[9]shows the range of light-curve rise 
times with viewing direction. As with the toys models, it can be 
seen that there is significant diversity in the rise time and that it 
tends to be shortest when viewing down n (i.e. close to ^ = 0) 
and longest when viewing closer to — n (i.e. 9 ~ n). In the next 
section, we examine the relationship between rise time and peak 
magnitude in some detail. 

5 IMPLICATIONS OF THE 3D CALCULATIONS 

We now consider some ramifications of the results obtained from 
the 3D model that has been discussed in Section|4] 



Figure 9. As for Figure [T] but showing the time at which maximum hght 
occurs. 
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Figure 10. The peak magnitude (Mp) versus the light curve rise time (tp) 
from the 3T2d200 model (Ropke et al. 2006). The opacity depends on com- 
position via equation 5. The three curves - dashed, solid and dotted lines - 
show the expected Arnett relation obtained from equations 3 and 4 adopt- 
ing a = 0.8, 1.0 and 1.2, respectively. The shaded grey-region indicates 
the ran ge of values for Mp predicted following .Stritzinger & Leibundguil 
)2005h for a ^^^Ni mass of 0.488 Mq. 



5.1 Amett's Rule and observational estimates of Ni masses 

In Section [X2I we used the toy models to investigate the possible 
systematic effects of an off-centre explosion on the inference of Ni 
masses. Here, we extend that discussion using the results obtained 
with the 3T2d200 explosion model. 

Figure [Tol shows the peak magnitude versus the time of max- 
imum light for the set of 100 light curves obtained for different 
viewing directions in the 3T2d200 model. As in Figure |4l curves 
are over-plotted to indicate the trends predicted by equations 3 and 
4. Similarly, the grey band indicates the region in which the com- 
put ed values of Mp are consistent with the relationship presented 
by jStritzinger & Leibundgutll200 5^. 

The figure reveals several interesting phenomena. Amongst 
the highest luminosity and fastest rise-time light curves there re- 
mains a trend which is reminiscent of what was suggested by 
the toy models (Figure However, this trend does not continue 
through the light curves with longer rise times; instead, the trend 
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first flattens and then turns over - a result of the non-hnear angular 
dependence of Mp as shown in Figure|7] This has the consequence 
that, as the viewing angle is varied, the light curve properties stay 
fairly close to those inferred via Amett's Rule (with a ~ 0.8) but 
with systematic variation of around ~ 0.2 mag (~ 20 per cent in 
luminosity). 

The value for the peak bolometric magnitude obtained us- 
ing th e model ^'^Ni mass and following Istritzinger & LeibundgutI 
( l2005l) is very close to the median of the peak magnitude ob- 
tained amongst the different viewing angles. The error range which 
they suggest - based on an estimated range of possible light curve 
rise times - happens to be rather close to encompassing the range 
of peak magnitudes obtained from the model. However, there are 
a moderate number of model light cur ves which lie outside this 
range, all of which are brighter than the lStritzinger & Leibundgu^ 
( 12005 ^ relation would anticipate. 

Thus we conclude that the calculations involving the fully 3D 
model support the primary result obtained with the toy models, 
namely that observationally significant systematic errors can be in- 
troduced via viewing angle effects in off-centre explosions. The 
particular model used here has a maximum spread in A/p of about 
0.5 mag, which is comparable to that obtained with toy model A. 
However, the particularly simple, linear dependence of A/p with 
cos 9 (see Figure|2} is not preserved in the 3D model. 

5.2 Scatter in the local Hubble diagram 

In this section we investigate what effect the predicted viewing an- 
gle dependence would have on the observed relationship between 
peak magnitude and recession velocity for SNe la in the Hubble 
flow. 

Assuming that the properties of the 3T2d200 model (Figure[7]l 
are representative for a population of SNe la, the viewing-angle de- 
pendence of the bolometric light curve would lead to a scatter about 
the mean relationship between peak flux and recession velocity. To 
quantify this effect, we have used the probability distribution of 
the peak magnitude (Figure[8} to map the region of the local Hub- 
ble diagram that would be occupied by a hypothetical population 
of SNe la which are in the Hubble flow and which have proper- 
ties as predicted by the 3T2d200 model. This region is indicated 
by the shaded area in Figure[TT] Over-plotted are lines which show 
contours of the probability distribution; these contours are slightly 
concentrated toward the fainter side owing to the asymmetry of the 
probability distribution (see Section l43] l. 

In constructing the Hubble diagram (via equation 10 of 
Stritzinger & Leibundgut 2005) we have adopted a value for the 
Hubble parame ter, Hp = 85 km s~^ Mpc~^. This value is the one 
that Stritzinger & LeibundgutI ( l2005h obtain from their analysis us- 
ing a typical ^®Ni mass of 0.42 M©; as they discuss, this value 
for Ho is unusually large and most likely a result of a net under- 
prediction for the amount of ^^Ni produced in the pure deflagration 
explosion models that they discuss. Although this is clearly an im- 
portant issue, it is not of relevance to our discussion of scatter in 
the Hubble diagram since the effect of changing the adopted value 
of Ho amounts to only a fixed offset in the ordinate of Figure [TT] 

For comparison, the flux and recession velocities for the sam- 
ple of twelve local SNe la which was assembled and used by 
Istritzinger & L eibundgut (2005) are also shown in the figure. Al- 
though relatively small compare d to data sets obtained from the 
large scale supernova surveys, the lStritzinger & Leibundgut ^20051) 
sample is ideally- suited for comparison with our results since the 
objects they include have been selected on the basis of having suf- 
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Figure 11. The peak UVOIR flux (F) versus recession velocity (cz) for 
SNe la in the local Universe. The grey shaded area indicates the region 
which would be populated by a multitude of 3T2d200 SNe in the Hubble 
flow having random orientations relative to the observer line-of-sight. The 
solid lines are drawn to delineate equal intervals of the underlying prob- 
ability distribution for Mp (they show 20 per cent intervals). The twelve 
data points which are over-plotted are those for the SNe la discussed by 
.Stritzinger & Leibundgut, ( ,2005iV 



ficiently large wavelength coverage that the bolometric peak mag- 
nitude is well-determined observationally; this is important since 
our calculations do not distinguish between different photomet- 
ric ba nds. Also, the fluxes tabulated by Stritzinger & Leibundgut 
(l2005h have not been normalised or modified using any empirically 
motivated stretch factors or light curve shape correlations. Since we 
do not compute band-limited light curves, we ca nnot use the st an- 
dard empirical correlations quantitatively (e.g. the |PhilliDsll993l re- 
lation) to renormalise our calculations. Therefore, for quantitative 
comparison of like with like, we require the distribution of observed 
bolometric fluxes corrected only for observational effects such as 
reddening. 

Figure[TT]clearly demonstrates two important points which we 
now discuss in turn. First, it shows that the viewing angle variation 
predicted by the 3T2d200 model does not substantially exceed the 
scale of scatter in the observed Hubble diagram. Therefore, one 
cannot rule out the 3T2d200 model as representative of a typical 
SN la by arguing that its asphericity is inconsistent with the ob- 
served narrowness of the Hubble-diagram correlation. This is in 
contrast to the more extreme toy models (see Section U!2l ( which 
could be excluded by this argument directly. 

Secondly, it can also be seen that the angular variations in 
the peak magnitude are on the correct scale to account for the 
observed dispersion about the mean relation. Thus it is proba- 
ble that orientation effects in aspherical SNe la could contribute 
significantly to the dispersion in the Hubble diagram. Indeed we 
conclude that, excluding the two apparently sub-luminous objects 
(SN1992bo and SN1993H), the width of the distribution of bolo- 
m etric magnitudes in the sample of well-observed objects selected 
by Istritzinger & LeibundgutI ^20050 can be explained fully by the 
viewing angle effects predicted from the 3T2d200 model without 
the need to invoke any other systematic or random uncertainty. 

If the scatter found in the bolometric model light curves trans- 
lates into a similar scatter in band limited light curves, this effect 
would have important implications for the use of SNe la as stan- 
dard candles: as a real physical effect which is dependent on the 
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explosion mechanism, it must be understood and cannot be elimi- 
nated as a source of dispersion by improved observational statistics 
or measurements alone. In addition, the finding that the probability 
distribution is not symmetric around the mean (cf. Figure[8l( could 
potentially introduce an additional bias in the interpretation of ob- 
served SN samples - if a significant fraction of the observed SN la 
sample has properties similar to the ones found in this model, the 
probability to find objects that are fainter than the mean may be 
somewhat larger than finding ones brighter than the mean. 

The commonly used methods to correct the peak magni- 
tude based on the shape of the Ught cu rve in certain filter bands 
(e.g., |Phillipslll993l iRiess et al.l I l996h are not expected to ac- 
count for the orientation effect fully. The calibration methods 
address a scatter in the band-limited peak flux which is most 
likely introduced by different degrees of nuclear burning and mix- 
ing which affect the emerging radiation in specific filter bands 
(Hoflichetal. 1995; Pinto & Eastman 2001; Mazzah et al. 2001; 
[Kasen & Wooslevl2007h through the resulting variation in the com- 
position and thermal structure of the ejecta. Because it is unlikely 
that the orientation effects follow the same correlation, it is ex- 
pected that these two sources of scatter will be superimposed in the 
distribution of observed peak magnitudes in particular photometric 
bands. 

It is conceivable that a relation between observables can be 
found that allows a correction following the spirit of the Phillips re- 
lation. However, the calculations presented here are not suitable to 
address this question because band-limited light curves, and likely 
also simulated spectra for the 3D model, are required to consider 
this point in detail. Knowing how the orientation effect affects the 
band limited light curves would also enable us to analyse the much 
larger data sets that are available for SNe la observed in specific 
wavelength bands to estimate the significance of asphericity in the 
observed objects. 



6 SUMMARY 

Motivated by recent explosion models involving off-centre ignition 
in SNe la, we have investigated the possible range of effects of a 
lop-sided distribution of nuclear ash on bolometric light curves for 
such objects. 

By studying results from both a grid of simple toy models and 
one real explosion model, we find that off-centre distributions of 
material which has undergone nuclear burning are likely to leave 
detectable imprints on observed light curves. An angular depen- 
dence of the light curve peak brightness is introduced - based on the 
models considered, the scale of this effect can be readily ~ 0.2 mag 
and conceivably much greater under extreme circumstances. 

This effect is large enough to have ramifications for the inter- 
pretation of the diversity in observed SN la properties and we have 
shown that the explosion model which we considered teopke et al.l 
|2003) already predicts sufficient viewing-angle sensitivity to ex- 
plain much of the typical scatter in a sam ple of objects with w ell- 
observed bolometric properties jStritzinger & Leibundgulll2005h . 

In view of the potential observational significance, further 
study of this class of explosion model is required. Investigation of 
viewing angle effects across a more diverse range of observable 
quantities is warranted (i.e. photometric band-limited light curves 
and, ultimately, spectra) to establish whether and by what means 
such effects may be observationally distinguished from other fac- 
tors which determine light curve properties (e.g. the total ^^Ni 
mass). Furthermore, a wider investigation of the possible extent and 



frequency of off-centre ignition in explosion models is needed so 
that the likelihood for occurrence of such events - and therefore the 
potentially important systematics they may introduce to statistical 
analysis of the SN la sample - can be quantified. 
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